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Abstract 

In this contribution, a stochastic theory for a branching process in a 
neutron population with two energy levels is investigated. In particu- 
lar, a variance to mean or Feynman-alpha formula is derived in this 
generalized scenario using the Kolmogorov forward or master equa- 
tion theory for the probabilities in a system with a compound Poisson 
source. 



1 Introduction 

There exist several relatively new applications where the energy distri- 
bution of the neutrons plays a significant role. One particular case is a 
method used in nuclear safeguards, namely the stochastic generalization 
of the so-called differential die-away analysis (DDAA) |T,'2','3','4]. Tradition- 
ally, the DDAA method was used as a deterministic method of detecting fis- 
sile material embedded in moderating surroundings using a pulsed source. 
The newly explored method, called differential die-away self-interrogation 
(DDSI) utilizes the inherent spontaneous neutron emission of the sample 
[5]. In the DDSI method the temporal decay of the correlations as a func- 
tion of the time delay between two detections of fast neutrons is used. 
This corresponds to a Rossi-alpha measurement with two energy groups. 
Furthermore, in recent pulsed experiments measuring the reactivity in fast 
cores of accelerator driven sub-critical systems it is found that two expo- 
nentials appear, indicating that the temporal behavior of the fast and ther- 
mal neutrons is separated in fast reflected cores. This leads to the fact 
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that a two group versions of the Feynman and Rossi-alpha formulas are 
needed [6]-[7]. 

In this contribution, a stochastic theory for a branching process in a 
neutron population with two energy levels is investigated based on the pre- 
vious results in Refs [6J- [7]. In particular a variance to mean or Feynman- 
alpha formula is derived in this generalized scenario using the master 
equation or Kolmogorov forward approach. The model includes absorb- 
tions, down scattering from fast to thermal neutrons, thermal fissions, de- 
tections and an external source of fast neutrons. Higher moments will also 
be discussed as well as specific applications to areas within safeguards 
research as well as specific applications to areas within safeguards re- 
search. 

2 The variance to mean via tlie forward Kol- 
mogorov approach 

In this section we will describe the two particle type system by using 
the Kolmogorov forward approach. We will include a compound Poisson 
source of fast neutrons described by the source strength Si which releases 
n particles with probability at an emission event (i.e. spontaneous fis- 
sion). The source is assumed to be switched on at time t = to, although 
dependence on to will not be denoted. The effects of detecting particles 
will also be included, denoted by the intensity A^. We will start by giving the 
differential equation for the probability P{Ni,N2,Zi,t) for Ni fast, N2 ther- 
mal neutrons at time t and Zi detected fast particles in the interval (0,t). 
We have summed all mutually exclusive events during an infinitesimally 
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small time interval dt and find for probabilities, 

+ Xia(N, + l)P(N, + l,N2,Z,,t) 
+ X2a{N2 + l)P{N,,N2 + l,Z,,t) 

+ XR{N, + l)P{N, + l,N2-hZ,,t) 

+ X2f{N2 + l)J2fik)P{Ni-k,N2 + l,Zi,t) 

k 

+ XdiNi + l)P{Ni + l,N2,Zi-l,t) 

+ 5i^p,(n)P(iVi-n,7V2,Zi,t). (1) 

n 

Here, Ai and A2 are the decay constants (total reaction intensities) for fast 
and thermal neutrons whereas Ai^, X2a are the absorbtion intensities of 
fast and thermal particles, respectively. The removal of fast particles into 
the thermal group is described by Xr while fission resulting from the ther- 
mal particles happens with the intensity of A2/. The intensities are related 
through 

Ai = Ai„ + Xr + Xd, (2) 

and 

A2 = A2a + A2/. (3) 

We will now solve this differential equation by using the generating function 
of the form 

G{X,Y,Z,t) = J]J]J]X^^y^^Z^^P(7V,,7V2,Zi,i), (4) 

Ni N2 Zi 

and describe the process in the time evolution of the generating function 
as, 

— = (X,a + XRY + XaZ-X,X)— + {X2a + X2MX)-X2Y) — 
+ S^{r{X)-l)G, (5) 
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where 

u{X) = Y.hX\ and (6) 

k 

r{X) = (7) 

n 

Here, fk is the probability of having exactly k neutrons produced in a fission 
event. Differentiation of equation ^ with respect to (X, Y, Z) and then 
letting [X = Y = Z = I) yields differential equations for the expectations 
as, 

^(iVi) = -Ai(iVi) + A2/z/i(iV2) + 5iri, (8) 
^(iV2) = -A2(iV2) + Afl(iVi), (9) 

= A,(iVi). (10) 

Here we have used the definition of the derivatives on the equations (HJ) 
and as z/i = dq/dX\x=i and n = dh/dX\x=i. We note that due to the 
source term with intensity Si the dynamical system consisting of equations 
^ - ^ will reach a steady state 0) and we find the stationary 

solution, 

A1A2 — ViARA2f UI1UI2 

(iV2) = iV2 = Ml!l, (12) 
(Zi) = eA2/iVit, (13) 
where e = Ad/A2/ and we have used the additional definitions ui and 002, 

-oji = -^(Ai + A2) + ^ V'lAi - Aa)^ + 4AiA2i/e//, (14) 

-W2 = -^(Ai + A2) - ^ V'^^^ ~ ^2^' + 4AiA2l/e//, (15) 

A1A2 



The expectation of the detections is found by integrating equation (1 0) and 
we note that the number of detections increases linearly with time. In order 



4 



to find the variance of the detector counts we need to determine the sec- 
ond factorial moment by yet another differentiation with respect to (X, Y, Z) 
and then letting {X = Y = Z = 1). The variance of the detector counts 
can be determined through the relation cr| = (Zi) + fizz where the mod- 
ified variance nzz is defined as nzz = {Z{Z - 1)) - (Z)^ = cr|^ - (Z) 
while in general we have fixv = {XY) - {X){Y). The differentiation pro- 
cedure gives a system of six dynamical equations of the modified second 
moments as 

d 

^/ixx = -2\ifixx + 2uiX2ffixY + i^2>^2fN2 + 5'ir2, (1 7) 

d 

^/ixy = -(Ai + A2)/ixy + A/j/ixx + '^iA2//iyy, (18) 
d 

it-Hyy = -2A2/iyy + 2AijyUxy, (19) 
at 

d 

g^f^zx = -Ai/Uzx + '^iA2//Uzy + AdyUxx, (20) 
d 

Q^f^ZY = -><2f^ZY + >^Rf^ZX + >^df^XY, (21) 

d 

Q^f^zz = 2eA2//Uxz, (22) 

where we have used the additional notations z/2 = d'^q/dX'^\x=i and r2 = 
x=i- Although, the system of equations (|1 7) - (|22|> is rather com- 
plicated and an analytical solution would be hard to find, we note that in 
stationary state the system breal^s down into two systems independent of 
each other where the moments jixx = fixx, ij-xy = fixY and /xyy = /2yy 
are constants. However, the equations describing detected particles need 
to be solved by e.g. Laplace transforms of ([20) and ^T} whereas it is pos- 
sible to find the sought moment fizz by integration by using equation (|22) . 
We find the constant 2nd modified moments as, 

(A^ + u;iu;2)('^2A2/iV'2 + Sir2) 

f^xx = ^TT — — , (^o) 

2(Ai + A2)wiW2 

A2AR(z/2A2/iV2 + S'ir2) 

f^XY = ^TTT — — r^; , (^V 

2(Ai + A2)c<;ia;2 

- A|(z/2A2/iV2 + Sir2) 

2(Ai + A2)UJiUJ2 
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The objective now is to solve (|20]) and ([21} by Laplace transform methods 
and we find the transformed identity as, 

sH[s) sH{s) 

with 

H{s) = s"^ + {u2 + uji)s + U1U2. (27) 

Note that we have assumed that the initial values of the moments fixz and 
fiYz were equal to zero at t = (at the start of the measurement), hence 
the roots of H{s) determine the temporal behavior of the Feynman-alpha 
formula. IVIoreover, the solution has many similarities to that found in Ref. 
(8]. The variance azz = (Z) + nzz'^s now found by integration of (|22]) and 
after some algebra the Feynman-alpha formula can now be written in the 
form 

Here, the complete expressions for Yi and Y2 are quite lengthy. However, 
it turns out that the sum Yg = + Y2 takes a rather simple form that also 
determines the value of the Feynman-alpha for large measurement times 

T 00 as, 

Yo = Y^ + Y2 = g2^^^^44^■ (29) 

We will now turn our attention to some quantitative examples of the Feynman- 
alpha formula in the form of Equation ([28)>. 



3 Results and discussion 

The Feynman-alpha formula for a two particle type system found by using 
the Kolmogorov forward approach including a Poisson source and effects 
of detecting particles is shown in Figure 1 (A and B). We have used the 
parameters ui = 3.0, U2 = 5.0, Si = 1.0, n = 1.0, r2 = 0.0, A2/ = 3/5 and 
Xd = 0.1. In Figure 1 A Ai^ = 2/3 whereas in Figure 1 B Ai = 1.0 and A2 = 
2.0. As expected, the curves grow according to Equation ((281) exponentially 
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to a maximum value determined by the constant Yq. However, this value 
can signifcantly vary depending on the intensities (Ai, A2, etc) involved the 
process. The increase of the curves is determined by two exponentials. 
Furthermore, unlike the case of the DDSI method, the presence of the 
two exponentials is not visible to the bare eye. In Figure 1A, it is seen 
that the ratio of the decay intensities for fast and thermal particles has a 
nontrivial effect on the maximum value by changing the ratio in the range 
(0.25-4.0). In Figure 1 B, the effect of the thermalization process described 
by the intensity Xr on the results is illustrated. Increasing thermalization 
increases the asymptotic value of the Feynman-alpha. 



Figure A Figure B 




t (arbitrary units) • (arbitrary units) 



Figure 1 : (A and B) The Feynman-alpha expression is shown for the pa- 
rameters z/i = 3.0, z/2 = 5.0, 5*1 = 1.0, ri = 1.0, = 0.0, A2/ = 3/5 and 
Xd = 0.1. In Figure A, Xr = 2/3 whereas in Figure B Ai = 1.0 and A2 = 2.0. 



4 Conclusions 

We have developed a forward Kolmogorov approach for the two group the- 
ory of the Feynman-alpha method, including a compound Poisson source 
and the detection process. The results agree with those calculated by the 
backward approach as reported in [6J. It is seen that, unlike in the DDSI 
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method (i.e. the two-group version of the Rossi-alpha method), the pres- 
ence of two exponents in the solution is not clearly visible. This means 
that detection of the presence of fissile material may not be as obvious 
as with the Rossi-alpha method. On the other hand, the determination of 
the exponents ui and u)2 by curve fitting could be more accurate in certain 
cases than with the DDSI method. However, the diagnostic value of the 
exponents in terms of determination of the sample parameters is not clear 
yet, and it requires further investigations, which will be reported in further 
worl<. 
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